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We derive an exact expression for the single particle Green function in the Calogero-Sutherland 
model for all the rational values of the coupling f3. The calculation is based on Jack polynomial 
techniques and the results are given in the thermodynamical limit. Two type of intermediate states 
contribute. The first one consists of a particle propagating out of the Fermi sea and the second one 
consists of a particle propagating in one direction, q particles in the opposite direction and p holes. 



I. INTRODUCTION 



The Calogero-Sutherland model describes interacting particles on a circle . Recent interest in the model arose 
from its relation with the Haldane-Shastry chain (/3 = 2) and also with the random matrix theory. Moreover, it has 
proven to be a good model to study fractional statistics. 

There has been several works devoted to derive the correlation functions in the Calogero-Sutherland model 
pT| , ^2| J^J^,p|,[7|,^0| . In particular, the density-density correlation functions were obtained. The Green function where 
the annihilation operator acts before the creation operator was also computed. Recently, Zirnbauer and Haldane ob- 
tained the Green function with the creation operator acting first for the special values (3 ~ 1/2, 1, 2 In this paper 
we compute this Green function in a simple way for all the rational values of (3 using Jack polynomial techniques. 

■ The method we use is closly related to our preceeding paper g, but the computation is more involved. There are 
t-H | also new physical phcnomcnas: Contrary to what one might naively expect, the intermediate states do not consist 

■ only of one particle which propagates out of the Fermi sea. There is an additional contribution where one particle 
' propagates in one direction, q particles in the opposite direction and p holes. 

, The model is defined on a circle of length L and the Hamiltonian given by : 

^ ' N 1 d 2 2 1 

D , where we use units in which % = m = 1. Its eigcnfunctions are known and written : 



*(a*) = ^(xMxi) (2) 

with : 



A(.)=n-( !Li ^) (3) 

and $(xj) is a symmetric function of the variables x\. The functions Q(xi) are symmetric polynomials in the variables 
Zj = e l2 ^ x i/ L anc i they are labeled by a partition A. They are known as the Jac k poly nomials J\{zi\ (3). The properties 
of these polynomials were recently studied in the mathematical literature [p^ , |18|Jlq 1 . 

The emergence of fractional statistics is seen by considering the spectrum of H . This spectrum is like a free spectrum 
and is described by a set of quasi- momenta ki . The momentum and energy are given by additive laws : 

9 N o 2 N 

Q a = ^J2 k » ^ = 9rE fc *> ( 4 ) 

i=l i=l 

and interaction between the particles is encoded in the fact that the quasi-momenta ki obey a generalized exclusion 
principle : 



H+l 



+ Ii. (5) 
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where Ii are positive integers or zero. The ground state is described by the configuration of fej's the most densely 
packed around the origin, I L = 0. For N odd, the ground state momenta are : 



fcf>=/^ N + 1 



By analogy with the fermion case ((3 = 1) we call this configuration the Fermi sea. To describe the elementary 
excitations it is convenient to multiply the quasi-momenta fcj by q so that they differ by integer numbers. Then, the 
Fermi sea can be described by a set of occupied quasi-momenta separated by p — 1 unoccupied ones. In this occupation 
number representation, a particle corresponds to one 1 followed by p — 1 zeroes and a hole corresponds to a sequence 
of q zeroes. Note that the effect of removing q particles from the Fermi sea is to create p holes (see figure I). 

a ) ... 0000(1 0)(1 0)(1 0)(1 0)(1 0)(1 0)(1 6) 000 00 00 0000.. . 

b) ... (1 0)0 0(1 0)(0 )(1 0)(0 )(1 0)(1 6) (1 0)(1 0) ... 

Fig.l:a) The ground state at j3 = 2/3 b) an excitation with 3 particles and 2 holes . 

This description does not tell which type of excitations are produced when one acts with a specific operator. This 
information can only be obtained by considering the properties of the wave functions. In the next part, we use these 
properties to obtain the Green function form factor. 

II. GREEN FUNCTION. 

The quantity of interest here is the Green function with the creation operator acting before the annihilation operator. 
The results in this section follow the paper ||] where the notations are defined. It is written : 

A r<0|*(x,t)* t (0,0)|0> A r = £, -<0\n^t)\ a > N + 1N + 1 < a \^(0fl)\0>, 



"> jv+i<a|a>jv + lN<0|0>jv 

jv<o|^(o,o)|q>jv+i| 2 

'\a> JV + i<a|a>jv + lN<0|0> JV ' 



V |w<0|^(0,0)|a>w + ir r -iE„t-iQ„x (a\ 

^\a> w + 1 <a|Q> w + i W <0|0> w e W 



where \a >n+i is a complete basis of states with N + 1 particles and |0 >n has N particles. This complete set of 
states can be expressed in terms of Jack polynomials in N + 1 variables : 

JV+l 

\a > N +i = \n, A >= Y[ ZiJ\(zi;P)\0 >n+i (7) 

i=l 



with 



N+l 

\0> N+1 =U^ {N+1)/2 U^-^- ( 8 ) 

i— 1 i<j 



and 



io> w =n^ (jv - i)/2 n^-^- o) 

i<j 



n is an arbitrary integer and the partition A is such that l(X) < N + 1 in order to avoid double counting of states. 
The quasi-momenta fcj entering the expression of the energy and momentum (4) of the state \a > are given in terms 
of A by : 
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ki = \ N - i+1 +n + p(i-^-^j . (10) 

The vacuum state |0 >n+i has been chosen in such a way that the quasimomentum of the N + 1 particle is at the 
left of the Fermi sea of the N particle vacuum. Another possible choice would have been to put this particle to the 
right by replacing the exponent of Zi by —f3(N — l)/2 in the expression (8) of |0 >at+i- If is an integer, these two 
states belong to the same Hilbert space and the two definitions give the same form factor in the thermodynamical 
limit. If (3 is not integer, we shall see that the two possible choices give different form factors in the thermodynamical 
limit. 

The action of ^(0,0) is to remove a particle at the positions Xi = 0, therefore setting the variable Zi to 1. The 
matrix element we need to compute is : 

M a = N <0\Y[z?- f3 (l - Zl fj x (l,z fl p)\0 >JV , (11) 

i 

for all A. It is possible to expand the Jack polynomial in N + 1 variables on polynomials in N variables : 

J A (1, Zi', P)=J2 J A/,(1; P)Ju{Zi] p) (12) 

V 

and here J\/ V (l; /3) is different from zero only when A covers v and when the difference between A and v is at most a 
box per column. When these conditions are satisfied, the coefficient is given by (is| : 

JX/vM = 3 Uw + D + aWJxlw + aW + l ); (13) 

The notations here are as follow: s = (i,j) is a box on a Young tableau identified by its coordinates 1 < j ' < A;. 
Given a box, s, a(s) is the number of boxes on its right and l(s) the number of boxes under it. The notation (...)a 
indicates that the quantities are evaluated with respect to the partition A and C'x/^ denotes the set of columns of A 
which have the same length as the corresponding column of v. 

The product in (11) can also be expanded on Jack polynomials with the relation ]l5| : 

Hil-Zif^^ifctiM^P) (!4) 

with coefficient H^(/3;fi) : 

g'(ftA0 = II i (15) 

f. A 1 + a(s) + pits) 



with notations defined above. The results of fl5|| are proven for f3 integer but it is simple to extend the proof to 
arbitrary rational /3 using Kaneko's generalized Selberg integrals ]l7j| . This coefficient is zero whenever the partition 
\x has more than p legs or more than (g — 1) arms for (3 = p/q. This is at the origin of the selection rule for the 
intermediate states. Using these expansions and the orthogonality relation for the Jack polynomials [jl6| : 

n < 0| J v {zi] (3) J P { Zl - /3)|0 >at= Su^pip, N) (16) 

we find the matrix element M a : 

M a = {-l) r3N YHP(f3;v + n)J x/v {l;P) - (17) 

where v + n is the partition v with n columns added (or subtracted when n negative) . In this expression v + n is a 
partition with at most p legs, therefore n < p. 
The correlation function is now written : 

N < 0\V(x, t)¥(0, 0)t |0 > N = Wc\ 2 ^ QaX ~ lEat - (18) 
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The evaluation of M a in the thermodynamical limit is our main result. All the quantities entering (17) are known 
from the mathematical literature and can be also found in ||. 

There is a first contribution coming from n — p and n = —r which will be denoted It this part, the only 

partitions contributing correspond to Young Tableaux, A, having one "arm" or one part, A = (Ai), with Ai an arbitrary 
positive integer or Young tableaux having N parts of length r denoted A = (r ). Each type of tableaux account for 
a particle propagating on a different side of the Fermi sea. In the thermodynamical limit, we obtain: 

/oo / -| \ 0— 1 

dw V^Ti) 2co <Q x ) e ~ iEt ( 19 ) 

where we used the variable : 

A. - — (20) 

E = ^^w 2 , Q = p /3™, (21) 

with pq — N/L. This is interpreted as the propagation of one particle over the Fermi sea. For p = 1, the fermions, 
we find the correct result. 

The contributions for n — 1, ...,{p — 1) are present for finite N but in the thermodynamic limit they are suppressed 
by factors of N. The only remaining contribution is for n < 0. It comes from partitions A which are of the form : 

(r + Ai, ...,r + X q , (r+pfp- q , (r + p - 1)^-^, ...,r N -^) (22) 



and 



J ^2 



Fig. 2: Young tableaux contributing to the correlation for (3 = 2/3. 

Here the A^'s {i — l,..,q) and r = —n are non negative integers and correspond to q particles propagating in the 
forward direction and one particle in the backward direction outside the Fermi sea . The X[ {i — 1, ...,p) are bounded 
by N and are related to the momenta of the holes inside the Fermi sea. 

When we take the continuum limit with the following variables : 

Xj- (/3- l)i _ Wj - 1 /3A- - i _ 1 - Vj r w + 1 Uj - (J3 - l)i _ & - 1 

f3N 2 ' (3N ~ 2 ' (iN ~ 2 ' j3N 2 1 ' 

we obtain the second part of the correlation function . The expression is given by : 

G^(x,t) = C(f3) J\dw 3 dw a Y[dVi Ffawj) e-u^-S-*, (24) 

Jl j=1 J -co J-i i= i 

The energy and momentum are : 

\ i=l i=0 ) \j=l i=0 / 
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and the constant : 

p 9-1 

2p\q\T{p/q) 

The form factor F(vi,Wj) is equal to ; 



K 



1W«*-»;>11 



(■Wi —Wj) 2 



\8-l 



Vl . ■•0 Vp 



n fc < 7 -(«*-«i)n?=i(nto( 



nSte-^)- 1 



(28) 



In the appendix we give the main lines of the derivation of IC and we argue that this expression can be simplified 
to a single contour integral : 



£ = (/3-l) 



r«(/?) 

2?ri 



dz ■ 



nLi(^-^) 



(29) 



where C w is a contour surrounding the points wi, w q as in the figure 3. 




Fig. 3: Contour of integration. The branch cuts of the integrand are represented by dashed lines. 

When (3 = 1/2, this expression of IC coincides with the result obtained by Zirnbauer and Haldane (T4J. For f3 integer 
(or, equivalently, q = 1) this result is rigorous and proves the conjecture in ]T^ |. The free fermion case j3 = 1 is special 
in the sense that the second contribution to the form factor vanishes. 

We tested numerically the equality of the integrals (28) and (29) for j3 = 1/2 and (3 = 3/2. 

This completes the computation of the form factor of the advanced Green function for all rational values of the 
coupling constant (3. 



III. CONCLUSIONS. 



We computed the advanced Green function in the Calogero-Sutherland model using Jack polynomial techniques. 
The computation done here follows from the techniques already used in 

The final expression for the form factor we obtain is the same for /3 integer or not. There is however a big difference 
between the two cases. If j3 is an integer, the form factor is invarient under a parity transformation ( change of sign 
of all the quasimomcnta which define the intermediate state a). If (3 is not an integer, one particle propagates in 
the forward direction and q particle propagate in the backward direction but there is no contribution coming from a 
particle propagating in the backward direction and q particles in the forward direction. The form factor is thus not 
invariant under a parity transformation. The origin of this phenomena can be traced back to two possible inequivalcnt 
iV + 1 particles vacua. Had we chosen the other possibility, we would have obtained the form factor conjugated under 
parity. The physical meaning of this phenomena is not clear to us. 

The results presented here simplify drastically in the continuum limit. Finding a way to do these computations 
directly in this limit would be of great interest. The complexity of the part described by K. in our results show 
that simple bosonic vertex operators don't seem to reproduce the right answer. Note however that the integrals over 
the variables £ in the first expression of IC bear striking similarities with the screening operators which occur in the 
Coulomb Gaz representation of conformal field theories 123] 
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APPENDIX A: SUMMATION. 



Let us go back to the expression for the form factor : 

M a = (-l)^£tf' s (/3;i/ + n)J A/ , l (l;j9) 



K(f3,N) 



^N a (f3;N)Nx(f3;N+l) 



(Al) 



and let us look at this expression for a fixed partition A and a negative integer n. Then the coefficient H^((3; v + n) 
forces the partition v + n to have p legs and (q — 1) arms. Since n is negative, the partition v must have r = —n 
columns and then p legs and (q — 1) arms. 

The main difficulty of the computation is to make the summation in (Al) on the partitions v such that J\/ v 7^ 0. 
As explained in Q , only cases for which A^ — v[ — 0, 1 are allowed. The case considered in figure 2 corresponds to 
(3 = 2/3. Again, we characterize the partition A by the numbers r, Ai, A^, which define the quasi-particle momenta in 
the continuum limit. In this case the partition A is of the form : 

(r + Ai, r + A 2 , r + A3, (r + 2) A ^ 3 , (r + l) x '^ ,r N ~^), 

and v can be read similarly from the figure : 

(r + V!,r + U2, (r + 2)^- 2 , (r + , r N ~^ ). 



Fig. 4: Young tableaux setting the notation. 

When we keep v\ fixed, we must sum over the possibilities v[ = X' t , A^ — 1 corresponding to the black boxes in the 
figure. The restriction that X/v is a one dimensional strip gives the following constraints on the i/j's : 



A2 < v\ < Ai, A3 < 2/2 < A2. 
If we use the variables appropriate for the continuum limit 

Xi - (P - l)i wi-1 pX'i-i l-Vi r w + 1 14 - ((3 - l)i & - 1 



(3N 



(3N 



r 



0N 



(A2) 



(A3) 



and we will use a = The sum over the iVs m the intervals defined by (A2) become integrals over the continuum 



variables £j's. This correspond to the shaded parts in the figure. 

Let us concentrate on the sum over i/, factorizing out the part which is independent of v[. This sum results in the 
following expression : 



G 



IN 



/ iei 



A W k -Vj-(P- T-r 1 £ k - V 

A A w k -Vi -1- J- 



fc=0 



n 



u,- - «i + a 



fc =i "/ 



(A4) 



Here J is the ensemble for which = v\ and we have to sum over all possibilities I C {1, ...,p}. The symbol |7| 
means the cardinal of the ensemble. Before taking the thermodynamical limit of this expression, we note the following 
property : 



iimly(-i)mrr 



f{vi - a) 

f(Vi) 



7)j — 



Vi + a 



Vi - Vi 



dv 1 -d Vp Ui<j( v i - v j) IlLi f(Vi) 



(A5) 



This expression, which comes from the definition of the derivative, can be also written as : 

uij...,i*p >\J IL\. ..lip J I V% I . 



/(«») 

/ iei J y 1 ieijei 



n 

ieljgl 



Vj -Vi+CLi 



Vi - Vi 



n 



Vi - Vi 



Vi - Vi 



Vi - Vi 



(A6) 



and it is unchanged if a quantity of order one in a, is added. Let f(v) have the form f(v) — J" 1 (v)/^ 2 (v). To first 
order in a; : 



f(vi-ai) h{vi- a 1 ai)f 2 {vi- a 2 ai) 



f(Vi) 



h(vi)f 2 {vi) 



= 1 - Oj [Oil 



f'M) , . /2(«i) 



a 2 - 



7l(«i) /2(«i) 



+ OK 2 )- 



Replacing this in (A5) and putting a\ = ... = a p = a gives the large N limit of the sum in (A3), with /i = IIiLi( w i — v), 
h = n?=ite - v ) and ai = 1 - /3, a 2 = -1. 



Q = 



d Vl •..0 Vp 






Uk<j( V k -Vj)U P j=l 







(A7) 



It remains to do the integrals over the variables £j's. When we multiply (A7) by the part depending explicitly on 
Vi and factor out the part which is independent of Vi we obtain the following multiple integral : 



di-Qp [n k<j (v k - v,) uu (nLoK - ntite - ^r 1 



(A8) 



Doing explicitly these integrals is not an easy task. However, a careful analysis of this expression suggests it can take 
the simpler form of a single contour integral. In the following we present the arguments leading to this conjecture. 

The integral I can be regarded as a complex function of the variables Wo, w q and v\, v p . We try to obtain the 
limit of I when the points wi, ...,w q collapse onto the point w and vi, ...,Vp — > v. Due to the antisymmetry of the 
integrand as a function of £i, changing counterclockwise Wi-\ and tUj, (i = 2, q) just multiplies the integral 
by a phase e n ^ 2fi ~ 1 \ We expect then branch points at Wi = wj. The singular part of the integral responsible of these 
branch points can be calculated: 



q-l - ■ 



dti q fl(£i-b) 9 flfl&-wi) f> - 1 

i<3 *=1 i=l 



r«cs) 



TV m II («>< - u>j) 2 " ^(wi, ...,«;,) 



(A9) 



The function R(wi, w q ) is symmetric, rational, homogeneous of order zero and it has no singularities, so it is 
constant. Evaluating the integral in the 'nested' limit w\ — ► (w 2 — ► ... — » (w 9 — » u>)...) gives R(wi, ...,w q ) = 1. 
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Further simplification arises when we specialize the variables Vi — v. This operation is delicate when done on an 
expression like (A7), so it is easier to restart from the discrete version (A4) and put Vi = v + ai. The sum reduces 
over the ensembles / = {1, ...s}, s < p : 



s=0 



n' Wk — v — ai — {(3 — l)a t-t — v — ai — a 
7/>7- — 7T — ni J- J- A. — 7! — m 



i=l Lfc=0 



Wk — V — ai 



fc=l 



35>S>'> 



s=0 



G(u + as) 



(A10) 



where 



k=0 



w k -v 



(I 



(3+1 



q-l 



n 



6 - u 



Here the notation (x)„ = r(x + n)/r(x) is used. In the last term of the equality (A10) one recognizes the discrete p th 
derivative of G(v). Using Stirling's formula and the homogeneity of G we get in the large N limit : 



(-lyp 



d" 



n?=o(^ 



II(« 

Li=0 



»=1 



(All) 



Note that the powers (3—1 and 1 have opposite sign in (A7) and (All). Now, we can take the limit of / when the 
Wi's collapse to w and all the Uj's are set equal to v. Without the constant and the factor n i<:j -(iOi — Wj) 2 @~ l , the 
result is : 



d p r 

(w - u; ) (/3_1)( ^ 1) (w - v )^ qP -^ [(«> " v) p -\w - vf/i- 1 



and this can be rewritten in the following form 



{(3 - l)(w - w) 2p ~ q (w - v)- pq (w - v)- p 



dP- 



dwP' 



(v - w) p 



(wq — w)P 



(A12) 



(A13) 



The equality of the two formulas can be checked by expanding the products under the derivatives and using the 
properties of the binomial coefficients. Then we write the derivative as a contour integral : 



dP- 



dwP~ 



(v — w) p 



(wq — w)! 3 



(P-1)! 
2m 



dz 



(v — z) T 



(z — w)p(wq — z)P 



The conjecture is that we can undo the specialization of the variables Wj and Uj to obtain 



1=03-1) 



2ni 



n K-^^nnK--^)- 1 

0<i<j i=l j=0 



dz 



IE=l(«i-*) 



(A14) 



(A15) 



where C w is a contour in the complex plane encircling w\, ...,w q . This function has the good analyticity properties 
we can expect from (A8). 

For (3 integer (or q = 1) this result is rigorous. In this case, there is no variable of integration £ and Q has simple 
poles at Wi — Vj, the residu being a polynomial of order 2q — 1 = 1 in each Vi. This justify the passage from (A13) to 
(A15). (3 = 1/2 is another case to test our conjecture. The contour of integration C w can be deformed to obtain C Wo . 
As (3 < 1 we have to subtract the divergent contribution of the contour at infinity, we thus obtain : 



1 2 1 r l 
4 JJ- Wi - v J_, 



dy 



v-y 



V w o-y wC^i -y)( w 2 -y) 



- i 



(A16) 



which is the result Zirnbauer and Haldane obtained using a supersymmetric calculation MA 
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